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ABSTRACT 

The spectral shapes of the contributions of different classes of unresolved gamma-ray emitters can 
provide insight into their relative contributions to the extragalactic gamma-ray background (EGB) and 
the natures of their spectra at GeV energies. We calculate the spectral shapes of the contributions to 
the EGB arising from BL Lacertae type objects (BL Lacs) and flat-spectrum radio quasars (FSRQs) 
assuming blazar spectra can be described as broken power laws. We fit the resulting total blazar 
spectral shape to the Fermi Large Area Telescope measurements of the EGB, finding that the best-fit 
shape reproduces well the shape of the Fermi EGB for various break scenarios. We conclude that 
a scenario in which the contribution of blazars is dominant cannot be excluded on spectral grounds 
alone, even if spectral breaks are shown to be common among Fermi blazars. We also find that while 
the observation of a featureless (within uncertainties) power-law EGB spectrum by Fermi does not 
necessarily imply a single class of contributing unresolved sources with featureless individual spectra, 
such an observation and the collective spectra of the separate contributing populations determine 
the ratios of their contributions. As such, a comparison with studies including blazar gamma-ray 
luminosity functions could have profound implications for the blazar contribution to the EGB, blazar 
evolution, and blazar gamma-ray spectra and emission. 

Subject headings: galaxies: active - gamma rays: diffuse background - gamma rays: galaxies 


1. INTRODUCTION 

The gamma-ray sky as currently observed by the Large 
Area Telescope (LAT) on board the Fermi Gamma- 
Ray Space Telescope consists of resolved point sources 
(such as active galactic nuclei (AGNs), pulsars, and star- 
forming galaxies), transient gamma- ray sources (e.g., 
gamma- ray bursts), and the diffuse gamma- ray radiation 
comprised of emission from the Galaxy and the extra- 
galactic gamma- ray background (EGB). The origins of 
the EGB are, as yet, unknown; however, it is expected 
that unresolved, extragalactic point sources provide a siz- 
able contribution to the EGB. 

In its first year of taking data, Fermi observed 1451 
resolved point sources with significance greater than 5a 
(First Fermf-LAT catalog (1FGL); Abdo et al. 2010b), 
of which 573 were associated with blazars (AGNs for 
which the jet is closely aligned with the observer’s line of 
sight (Blandford & Konigl 1979)). Thus, blazars con- 
stitute the largest class of astrophysical objects asso- 
ciated with gamma-ray sources, and it has long been 
suspected that unresolved blazars should provide a sub- 
stantial contribution to the EGB (Padovani et al. 1993; 
Stecker et al. 1993; Salamon & Stecker 1994; Chiang 
et al. 1995: Stecker & Salamon 1996; Kazanas & Perl- 
man 1997; Chiang & Mukherjee 1998; Sreekumar et al. 
1998; Mukherjee & Chiang 1999; Miicke & Pohl 2000; 
Giommi et al. 2006; Narumoto & Totani 2006: Dermer 
2007; Kneiske & Mannheim 2008; Inoue et al. 2008; In- 
oue et al. 2010; Inoue & Totani 2009; Abdo et al. 2010c; 
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Stecker & Venters 2011). 

It should then, perhaps, come as no surprise that 
the Fermi-LAT, with its improved sensitivity enabling 
it to observe many more blazars, measured the inte- 
grated intensity of the EGB to be substantially lower 
than that measured by the Energetic Gamma-Ray Ex- 
periment Telescope (EGRET) aboard the Compton 
Gamma-ray Observatory in the 1990s (1.03 (± 0.17) x 
10~ 5 cm -2 s _1 sr _1 for Fermi versus 1.45 (± 0.05) x 
10~ 5 cm -2 s _1 sr -1 for EGRET; Sreekumar et al. 1998; 
Abdo et al. 2010g). However, exactly how much of the 
change in the EGB integrated intensity is due to the 
ability to resolve many more blazars remains unclear. 
On the one hand, since the Sreekumar et al. (1998) de- 
termination of the EGRET EGB, models of the galac- 
tic foreground emission have been updated to reflect re- 
cent observations of the interstellar medium (Abdo et al. 
2010g), accounting for at least some of the change. On 
the other hand, the determination of the distribution 
of blazars with respect to luminosity and redshift (the 
blazar gamma-ray luminosity function, GLF) and by ex- 
tension, their contribution to the EGB directly from their 
observed flux distribution (such as that performed in 
Abdo et al. 2010c) is non-trivial, for four reasons. 

First, blazars are variable at gamma-ray energies and 
without knowledge of the degree to which the gamma- 
ray luminosity of a blazar changes during the flaring 
state and the blazar duty cycle, the reconciliation of the 
observed blazar source counts with model GLFs is im- 
possible (Stecker & Salamon 1996). This is because the 
increase in flux of a blazar (sometimes up to an order 
of magnitude; see e.g., McLaughlin et al. 1996: Mukher- 
jee et al. 1997; Nolan et al. 2003: Vercellone et al. 2004) 
during the flaring period introduces a selection effect in 
that flaring blazars are easier to detect than quiescent 
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blazars, and quiescent blazars would have to be relatively 
bright and nearby to be detectable. The largest impact 
on the observed blazar source counts will be at fluxes just 
above the Fermi-LAT sensitivity as some blazars below 
the threshold make it into the sample because they are 
flaring. While the magnitude of the impact on the ob- 
served blazar source counts depends on the blazar GLF 
and the blazar duty cycle, the effect will be to flatten the 
faint-end slope of the observed blazar source counts (as 
blazars move from lower flux bins to higher ones) lead- 
ing to an underestimation of the blazar contribution to 
the EGB as extrapolated from the observed source counts 
alone. Abdo et al. (2010c) concluded that since the peak- 
to-mean flux ratio for most Fermi sources is ~ a factor 
of two, that there are no large systematic uncertainties 
due to blazar variability. However, a small peak-to-mean 
ratio is expected for blazars that are observed mostly in 
a flaring state, and as most of the observed blazars are 
fainter blazars that will be more represented by flaring 
blazars, then while it will be true that most blazars will 
have a small peak-to-mean ratio, it is not necessarily the 
case that the effect of blazar variability will be small. 

Second, the large angular resolution of the Fermz-LAT 
at lower energies (~ 5° at 100 MeV) could introduce sig- 
nificant source confusion for GLF models that predict 
high blazar densities resulting in fewer resolved blazars 
and a higher contribution to the EGB at lower energies 
(Stecker 1999; Stecker & Venters 2011). The impact of 
source confusion on observed source counts depends on 
the source density as given by the blazar GLF. As such, 
accounting for source confusion requires a priori knowl- 
edge of the source density, which is exactly what the 
observer wishes to determine. Thus, source confusion 
further complicates analyses of the blazar contribution 
to the EGB based solely on observed source counts (for 
a detailed discussion of source confusion and the differ- 
ence between the definition employed in this discussion 
and that employed by Abdo et al. 2010c, see Stecker & 
Venters 2011). 

Third, as demonstrated in Abdo et al. (2010c), the re- 
constructed fluxes of the individual sources are subject to 
a considerable amount of uncertainty, particularly at the 
faint end (see their Figure 6). Finally, the model of the 
Ferrm-LAT detection efficiency employed by Abdo et al. 
(2010c) assumed power- law spectra for the sources 0 . As 
such, these uncertainties complicate the determination of 
the distribution of blazars with respect to luminosity and 
redshift (the GLF) and the determination of the blazar 
contribution to the EGB from their observed flux distri- 
bution. Notably, Stecker & Venters (2011) found that 
contrary to the conclusion of Abdo et al. (2010c), the 
observed flux distribution of Fermi blazars does not, as 
yet., rule out the possibility that the EGB is dominated 
by emission from unresolved blazars (see also Abazajian 
et al. 2011). Thus, the contribution of still unresolved 
blazars to the EGB remains in dispute. 

Additional information about the contributions to the 
EGB can be obtained through studying the shape of its 
energy spectrum, and the shapes of the collective inten- 

5 It should also be noted that the model for Fermi-LAT detec- 

tion efficiency (particularly at energies below ~ 300 MeV where 
multiple scattering in the detector becomes important) is subject 
to change during the course of observations by Fermi-LAT. 


sity spectra of suspected contributors (Pavlidou & Ven- 
ters 2008, hereafter PV08; Venters et al. 2009; Venters 
2010). Analysis of the first year of Fermi data yielded 
an EGB spectrum consistent with a featureless power 
law with spectral index T ~ 2.4 (Abdo et al. 2010g). 
Upon first reflection, the Fermi EGB spectrum appears 
to be consistent with the hypothesis that the EGB is 
dominated by emission from unresolved blazars since the 
mean spectral index for blazars, Fq, is also ~ 2.4. How- 
ever, several effects complicate this simple picture. 

First, the astrophysical population of blazars is ac- 
tually composed of two separate sub-populations (flat- 
spectrum radio quasars, FSRQs, and BL Lacs) with 
distinct spectral properties (Fo ~ 2.45 for FSRQs and 
To ~ 2.2 for BL Lacs; Abdo et al. 2010c). Second, even 
within a given sub-population of blazars, the spectral in- 
dices of individual blazars form a distribution with some 
spread (Stecker & Salamon 1996; Venters & Pavlidou 
2007), which causes the collective spectrum of unresolved 
blazars to curve as harder blazar spectra become more 
important at higher energies (Stecker & Salamon 1996; 
Pavlidou & Venters 2008). On the other hand, this effect 
is somewhat mitigated by the spectral bias introduced 
by the fact that blazars with harder spectral indices are 
easier to observe in a flux-limited survey 5 6 , and thus, are 
more likely to be resolvable and not play as big a role in 
producing the EGB. Third, just as unresolved blazars are 
expected to contribute to the EGB, unresolved members 
of other known astrophysical gamma-ray emitters (such 
as star-forming and starburst galaxies) should also con- 
tribute, but with spectra that are substantially different 
from those of blazars (Pavlidou & Fields 2002; Fields 
et al. 2010; Lacki et al. 2011; Makiya et al. 2011; Stecker 
& Venters 2011) and may not even resemble power laws 
at gamma-ray energies. Finally, recent observations con- 
ducted by Fermi indicate that even blazar spectra can 
break at ~ GeV energies and would no longer be de- 
scribable by simple power laws (Abdo et al. 2009a; Inoue 
et al. 2008; Inoue et al. 2010; Inoue & Totani 2009). 

The intuitive conclusion would be that the combina- 
tion of these effects should cause the energy spectrum of 
the EGB intensity to exhibit features (Pavlidou & Fields 
2002). The observation instead of an EGB with a single, 
featureless power-law spectrum begs the question: is the 
lack of observed spectral features necessarily indicative of 
a single dominant source population with individual un- 
broken power-law spectra in the energy range of > 300 
MeV? If not, then what can such a featureless power- 
law EGB spectrum tell us about the relative contributions 
of the separate populations and how do they depend on 
their collective spectra? The answer to these questions 
provide novel constraints on blazar GLFs and thus could 
have profound implications for the cosmological prop- 
erties of known astrophysical gamma-ray emitters, their 
corresponding contributions to the EGB, and the general 
properties of their spectra at gamma-ray energies. 

In this paper, we calculate the collective spectrum 

6 It should he noted that the Fermi survey is not exactly a 
flux-limited survey due to the non-uniformity of the total diffuse 
background throughout the sky. As such, spectral bias is more sig- 
nificant in the Fermi survey than typical for a truly flux-limited 
survey. However, as evidenced from the analysis of the mean spec- 
tral indices performed by Abdo et al. (2010c), this appears to have 
more of an effect on BL Lacs than FSRQs. 
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of unresolved blazars with individual spectra exhibiting 
broken power laws. In so doing, we seek to determine 
whether such a population of unresolved blazars can re- 
sult in a featureless power-law spectrum resembling that 
of the EGB. If so, we could then investigate the impli- 
cations for the relative contributions of FSRQs and BL 
Lacs, which could, in turn, have implications for the cos- 
mological properties of blazars. In Section 2, we present 
the formalism of the calculation of the spectral shape of 
the collective unresolved blazar emission. In Section 3, 
we discuss the inputs of the calculation and their un- 
certainties. In Section 4, we present the results of the 
calculation, and we discuss them in Section 5. 


under consideration. For the first year of Fermi data, 
E m in ~ 2 x 10“ 9 photons cm -2 s~ J . 

Equation (4) can be characterized in terms of factors 
that determine the overall magnitude of the unresolved 
blazar contribution to the EGB (flux terms) and factors 
that determine the overall shape of the blazar contribu- 
tion (spectral index terms). For a carefully chosen defi- 
nition of p(F) (determined from the analysis of spectral 
indices of the flux-limited sample of Fermi blazars that 
accounts for the spectral bias inherent in a flux-limited 
catalog; see Section 3 and Venters et al. 2009), the mag- 
nitude and shape terms decouple, and Equation (4) can 
be rewritten as 


2. FORMALISM 

To calculate the collective spectrum of unresolved 
blazars, we follow the formalism as outlined in PV08 
with one major difference. Instead of taking blazar 
spectra to be simple power laws over gamma-ray ener- 
gies (F oc E~ r , where V is the photon spectral index 
at gamma- ray energies), blazar spectra are taken to be 
smoothly broken power laws: 


1(E) = 
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( 5 ) 

where 7o is a normalization constant depending on the 
flux distribution of unresolved blazars and 
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where Fe(E) is the differential photon flux in units of 
photons per unit area per unit energy per unit time, E b 
is the break energy, Ti is the low-energy slope, r 2 is 
the high-energy slope, and n quantifies the sharpness of 
the transition from the low-energy power law to the high- 
energy power law. For the purposes of this paper, we take 
n = 1. The total flux, F, of photons with energies greater 
than some fiducial energy, E;, is found by integrating 
Fe(E) over energy, 


F = F 0 
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For the purposes of this paper, we take Ef = 100 MeV. 
Then, the contribution of a single unresolved blazar to 
the EGB is 

F[(E/E b f^ + (E/E b f* n }- l/n 

1 1 ” r \oj 

4 tt /“ {(E/E b ) + (E/E b ) r * n ] /n dE 

where the flux of one source is uniformly distributed 
over the entire sky in anticipation of an isotropically dis- 
tributed cosmological population and I has units of pho- 
tons per unit area per unit energy per unit time per unit 
solid angle. 

Following PV08, we characterize the flux distribution 
of unresolved blazars as a function g(F) and the distri- 
bution of blazar spectral indices (or spectral index dis- 
tribution, SID) as a function p(F), where Fi = F — AFi, 
F 2 = F + AF 2 , and p( F), AFi, and AF 2 are determined 
from observations (see Section 3). Then, the total con- 
tribution of unresolved blazars to the EGB is given by 
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where F ni ; n is the sensitivity of the gamma-ray telescope 
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3. INPUTS 

As demonstrated in PV08, the unresolved blazar con- 
tribution to the EGB is not just a question of magni- 
tude, but also of spectral shape, and the spectral shape 
is sensitive to the distribution of blazar spectral indices 
at GeV energies. Both PV08 and Abdo et al. (2010c) 
assumed that blazar spectra at gamma-ray energies are 
power laws and did not account for possible breaks in the 
spectra. In this paper, blazar spectra take the forms of 
smoothly-broken power laws (see Section 2). 

For each sub-population of blazars, we determine the 
SID, p(F), from the likelihood analysis of Venters & 
Pavlidou (2007) fitting blazars from the Fermz-LAT First 
Catalog of AGNs (Abdo et al. 2010d) to Gaussian SIDs 
accounting for errors in measurement of individual spec- 
tral indices. Due to the survey bias towards harder 
spectral indices present in Fermi data, we applied the 
likelihood analysis only to the subset of blazars with 
photon fluxes > 7 x 10~ 8 photons cm -2 s _1 and galac- 
tic latitudes > 10° (as per Abdo et al. 2010c). We 
determined that the maximum-likelihood Gaussian SID 
can be characterized by a mean (Fo) and a spread (<to) 
with maximum-likelihood parameters determined to be 
Fo = 2.45 and <x 0 = 0.16 for FSRQs and Fq = 2.17 and 
<t 0 = 0.23 for BL Lacs, which are similar to the find- 
ings of Abdo et al. (2010c). Based on 1FGL spectra of 
two prominent blazars from the Fermi- LAP Bright AGN 
Sample (Abdo et al. 2009a, ; see Figure 1), we model the 
spectral breaks by taking ATi = 0.1 and AF 2 = 0 . 9 . As 
in the measured blazar spectra, we treat the break en- 
ergies distinctly for FSRQs and BL Lacs. We took the 
break energy, E b , to be 4 GeV for FSRQs and 15 GeV 
for BL Lacs. In so doing, we consider two cases: 

Scenario 1. - Blazars within a population evolve 

such that their break energies are observed to be 
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roughly the same. 

Scenario 2. - The break energies of the intrinsic 
spectra of blazars within a subpopulation are the 
same, but because of redshift effects, they are ob- 
served to be different. 

In order to perform the calculation in the Scenario 2, 
we separately model the contribution arising from dif- 
ferent redshift bins from an assumed GLF model for a 
given subpopulation of blazars and then determine the 
composite spectrum in the final step. For the purposes 
of this analysis, we assume the best-fit pure luminosity 
evolution model of Narumoto & Totani (2006) for FSRQs 
and their best-fit luminosity-dependent density evolution 
model for BL Lacs 7 . In this manner, we seek to inves- 
tigate the impact of changing the GLF model(s) of the 
blazars. We should also note that alternative statistical 
treatments on more extensive data sets as in Abdo et al. 
(2010f) could result in different estimates of the break 
energies and changes in spectral index. 

Another alternative treatment performed by Abdo 
et al. (2010c) consists of stacking the spectra of ob- 
served blazars. However, the analysis was performed 
on the flux-limited sample of Fermi blazars, which con- 
sists almost entirely of FSRQs (hence, the similarity be- 
tween the stacked spectrum for the flux-limited sample 
of blazars and that of FSRQs). Moreover, while such an 
analysis is effective in determining the average spectral 
index of a given population of blazars, it reveals nothing 
about the spectral properties of the population of blazars 
(e.g., the spread in the SID, spectral breaks). Further- 
more, as the Fermi-LAT survey is biased against sources 
with softer spectra, it is certainly biased against sources 
exhibiting significant breaks in their spectra. We also 
note that despite finding that the SIDs of blazars have 
non-zero spreads, Abdo et al. (2010c) did not account 
for these spreads in their calculated spectra of the blazar 
contributions to the EGB as evidenced by the fact that 
these spectra, as well as their estimated uncertainties, 
are still power laws (see Figures 18, 19, and 20 of Abdo 
et al. 2010c). Similarly, Abdo et al. (2010c) did not ac- 
count for spectral breaks in their calculated spectra of 
the blazar contributions to the EGB. Inclusion of these 
effects would introduce curvature in the spectra of the 
blazar contributions to the EGB and not result in the 
power laws indicated by Abdo et al. (2010c). 

As in PV08, we do not include information about the 
magnitude of the unresolved blazar contribution to the 
EGB. Instead, we normalize the collective spectrum of 
each sub-population in order to best fit (as determined 
using a x 2 -analvsis) the measured spectrum of the EGB 
from the Fermi- LAT first year data (Abdo et al. 2010g). 
For the composite spectrum of FSRQs and BL Lacs, we 
normalize the collective spectrum of each sub-population 
such that their total composite spectrum is the best-fit 
spectrum to the measured spectrum. 

We should note that while applying the aforementioned 
cuts to the sample accounts for much (though likely not 
all; see Venters et al. 2009) of the survey spectral bias, do- 
ing so also significantly reduces the sample size and will 

7 In so doing, we account for the observation that BL Lacs are 
on average situated at lower redshifts than FSRQs. 
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FIG. 1. — Sample broken power-law spectra for blazars. Top: 
PKS0235+164 measured spectrum with broken power- law (solid) 
and single power- law spectra (dashed; F = 2.1433). Bottom: 
3C454.3 measured spectrum with broken power-law and single 
power-law spectra (F = 2.4662). 


ultimately introduce more uncertainty in the determina- 
tion of the likelihood parameters. Furthermore, while in 
the flux-limited sample of blazars, spectral index and flux 
do appear to decouple (see Figure 1 of Abdo et al. 2010c), 
we expect it to be only approximately correct. A more 
detailed calculation, however, requires a statistical anal- 
ysis to determine the intrinsic blazar SID, which would, 
in turn, require knowledge of the blazar GLFs since the 
survey spectral bias depends on the redshift and lumi- 
nosity distributions of blazars® (Venters et al. 2009). In 
the case of FSRQs, applying the high-flux cut changes 
the SID very little (Abdo et al. 2010c; Stecker & Venters 
2011). so we do not expect the collective spectrum of un- 
resolved FSRQs to be appreciably different from that de- 
termined from the resolved FSRQs. On the other hand, 

8 For this same reason, we also neglect source confusion at low 
energies arising from the increase in Fermt-LAT point, -spread func- 
tion below 1 GeV. 
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in the case of BL Lacs, the sparseness of the population 
at fluxes above the cut introduces a considerable margin 
of error. As such, the collective spectrum of unresolved 
BL Lacs could be harder or softer than that expected 
from the resolved BL Lacs. As BL Lacs are, on aver- 
age, harder than FSRQs, they will likely contribute most 
significantly at the higher energies, and hence, the dif- 
ferences in spectra between resolved and unresolved BL 
Lacs will mostly impact the collective unresolved blazar 
spectrum at higher energies. 

Finally, we should note that recent multi-wavelength 
observations conducted by Fermi and other telescopes 
suggest that blazars can be further subdivided into the 
categories of low synchrotron peak, intermediate syn- 
chrotron peak, and high synchrotron peak (Abdo et al. 
2010e). In principle, one could gain further insight by 
applying this procedure to each subdivision as there is 
some indication that these categories are also spectrally 
distinct at gamma-ray energies (Abdo et al. 2010e). No- 
tably, as demonstrated in Abdo et al. (2010f), most high- 
synchrotron peaked BL Lacs (HBLs) do not exhibit spec- 
tral breaks in Fermi - LAT energy range while those HBL 
breaks that have been observed are different in character 
(going from soft to hard rather than hard to soft). How- 
ever, in light of the cuts on galactic latitude and flux, 
further subdividing the sample of blazars could lead to 
small sample sizes in some of the categories, and conse- 
quently to poor SID parameter determination. Thus, for 
the purposes of this analysis, we retain the original clas- 
sification scheme of FSRQs and BL Lacs. In any case, 
given that BL Lacs likely contribute most significantly 
at higher energies, the effect of the HBL spectra is likely 
to improve the fit at high energies, though a more de- 
tailed study is in order as more data become available. 
One might also be concerned about systematic changes 
in blazar spectral indices while flaring; however, analy- 
ses of EGRET blazar spectral indices found no evidence 
of such systematic changes in spectral index with flar- 
ing (Nandikotkur et al. 2007; Venters & Pavlidou 2007), 
and Fermi observations of individual blazars have thus 
far revealed no systematic changes in spectral index with 
time or flux (Abdo et al. 2009b, 2010a; Ackermann et al. 
2010). 

4. RESULTS 

The spectral shapes of the contributions of unresolved 
blazars to the EGB as determined from the Gaussian 
SIDs and broken power-law gamma-ray spectra discussed 
in Section 3 are plotted in Figures 2-4. In Figure 2, we 
have plotted the best-fit spectral shapes of the collec- 
tive intensities of unresolved FSRQs (x? e d = 4.1) and 
BL Lacs (x? ed = 3.3). In Figure 3, we have plotted the 
best-fit spectral shape of the total collective spectrum of 
unresolved blazars in Scenario 1 (solid line; x? ? d = 0.7: 
for explanation, see Section 3) and the individual con- 
tributions from FSRQs (dashed line) and BL Lacs (dot- 
dashed line). In Figure 4, we have plotted the best-fit, 
spectral shape of the total collective spectrum of unre- 
solved blazars in Scenario 2 (,\f ed = 0.6). For compar- 
ison, we have also plotted the best-fit spectral shape of 
the total collective spectrum assuming that FSRQ spec- 
tra break at Ei ~ 1 GeV and BL Lac spectra break at 
Et, ~ 5 GeV (assuming Scenario 1; Figure 5). Note that 
such spectra do not adequately fit the spectra in Figure 




Fig. 2. — Best-fit shapes of the individual collective intensities of 
the sub-populations of unresolved blazars. Top: the best-fit shape 
for FSRQs. Bottom: the best-fit shape for BL Lacs. Data points: 
spectrum of the EGB as measured by the Fermi - LAT (Abdo et al. 
2010g). 


1. As such, this alternative is included merely for the 
sake of comparison. 

As noted in Section 3, BL Lacs tend to be harder than 
FSRQs, and thus, the BL Lac collective spectrum is no- 
ticeably harder than that of the FSRQs. Also apparent 
is the effect of the higher break energy for the BL Lacs 
as compared with the FSRQs. As can be seen in Figure 

2, the collective spectrum of neither sub-population of 
blazars reproduces well the spectrum of the EGB: the BL 
Lacs are too hard while the FSRQs are too soft. How- 
ever. if the spectra are added together as in Figures 3 
and 4 (renormalized as described in Section 3), the re- 
sulting spectrum (solid line) reproduces well that of the 
EGB (being within the statistical error bars of nearly all 
of the data points and within ~ l.Scr of the last data 
point). The harder spectra and higher break energies for 
the BL Lacs compensate for the suppressed intensity of 
FSRQs at higher energies. Notably, the redshift depen- 
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FlG. 3. Best-fit shape of the total collective intensity of unre- 
solved blazars in Scenario 1 as described in Section 3. Dashed line: 
the contribution from FSRQs. Dot-dashed line: the contribution 
from BL Lacs. Solid line: the shape of the combined population of 
FSRQs and BL Lacs. Data points: same as in Figure 2 
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Fig. 4. — Same as in Figure 3, but for Scenario 2 as described in 
Section 3. 


dence of Scenario 2 shifts the breaks to lower energies 
resulting in the transition from FSRQ dominance to BL 
Lac dominance occurring at lower energies (~ 1 GeV 
rather than ~ 3 GeV). Thus, the relative contribution 
from BL Lacs to the total blazar collective spectrum is 
greater in Scenario 2 than in Scenario 1. In the case of the 
even lower break energies as shown in Figure 5, the total 
collective spectrum fits the EGB fairly well (\f cd = 1.1). 
but the relative contributions of the FSRQs and BL Lacs 
are such that the contribution from BL Lacs dominates 
at all energies. This is in contrast with the scenarios 
presented in Figures 3 and 4 in which BL Lacs dominate 
only energies greater than ~ few GeV. Also, the sharper 
contrast between Figures 3 and o than that between Fig- 



Fig. 5. — Same as in Figure 3, but for lower break energies. 


ures 3 and 4 indicate that the transition and the resulting 
relative contributions of FSRQs and BL Lacs are more 
sensitive to the break energies than model of the blazar 
GLF(s). This is due to the fact that each contribution 
is dominated by its closest and brightest (though still 
unresolved) members 'for which the observed break ener- 
gies are similar to their intrinsic values for the GLF(s) 
considered. 

5. DISCUSSION AND CONCLUSION 

We have calculated the spectral shapes of the contri- 
butions to the EGB arising from BL Lacs and FSRQs as- 
suming blazar spectra can be described as broken power 
laws. We found that in the case that blazar spectral 
breaks are indeed common neither sub-population alone 
can adequately reproduce the spectrum of the EGB. 
However, in a combined spectrum, the harder spectra of 
the BL Lacs could compensate for the softer spectra of 
the FSRQs, resulting in a collective blazar spectrum that 
is similar to that of the EGB within uncertainties. Fur- 
thermore, we have found that the relative contributions 
of FSRQs and BL Lacs required to fit the EGB spectrum 
are sensitive to the nature of the spectral breaks and, to 
a lesser extent, the blazar GLFs. 

The question of whether the relative contributions nec- 
essary to reproduce the overall spectrum of the EGB as 
determined in this method are reasonable given future 
Fermi measurements of the blazar GLFs (which account 
for the notable observational uncertainties discussed be- 
low and in Section 1) would provide insight into the over- 
all blazar contribution to the EGB as well as the nature 
of blazar variability. According to current Fermi obser- 
vations, the luminosities of FSRQs are greater than those 
of BL Lacs (L 7 ,fsrq ~ 10 2 x Abdo et al. 2009a, 

2010d). On the other hand, BL Lacs are likely situated at 
lower, on average, redshifts than FSRQs (Dermer 2007). 
Thus, it is possible that the closer proximity of BL Lacs 
(or their numbers) compensates for their deficit in lumi- 
nosity with respect to FSRQs resulting in their roughly 
comparable or slightly enhanced relative contribution to 
the collective blazar spectrum seen in Figures 3, 4, and 
5. However, it should be noted that the closer proximity 
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of BL Lacs and their harder spectral indices would also 
make them more easily observable by Fermi , limiting the 
contribution of unresolved BL Lacs to the EGB-though 
if the luminosity function for BL Lacs is broad enough, 
there could be many low-luminosity BL Lacs that would 
escape detection. 

It thus remains to be seen how the collective intensi- 
ties of the sub-populations of blazars actually compare 
with one another. In order for a true comparison to 
be drawn, the GLF for each sub-population of blazars 
needs to be measured. If it should be the case that the 
relative contributions are not reasonable given gamma- 
ray observations, then either the breaks included in the 
best-fit spectrum are not typical of blazars or blazars 
are not sufficient to explain the EGB. Already, Fermi 
studies of observed blazar flux distributions have sug- 
gested that emission from unresolved blazars may not 
comprise the dominant contribution to the EGB (Abdo 
et al. 2010c). However, given that the physics behind 
the gamma-ray emission of other known and speculated 
contributors (i.e., normal galaxies', cascades of ultra- high 
energy cosmic rays and TeV gamma rays, and dark mat- 
ter annihilation) result in spectra that are quite distinct 
from that of the EGB as measured by Fermi (see e.g., 
Ando et al. 2007; Kalashev et al. 2009; Siegal-Gaskins 
& Pavlidou 2009; Ahlers et al. 2010; Fields et al. 2010; 
Venters 2010; Berezinsky et al. 2011; Stecker & Venters 
2011), a close resemblance of the collective spectrum of 
unresolved blazars could be striking. 

The necessary reconciliation of the clues to the blazar 
contribution to the EGB provided by studies of the col- 
lective blazar spectrum with the clues provided by stud- 
ies of the blazar source counts could thus provide insight 
in the gamma-ray emission properties of blazars. Blazar 
variability plays a substantial role as flaring blazars 
would be more easily observed by Fermi than quiescent 
blazars. As such, studies of the blazar GLF, observed 
blazar flux distributions, and ultimately, the blazar con- 
tribution to the EGB are largely a question of the blazar 
duty cycle. If blazars spend the majority of their lifetimes 
in the more observationally challenging quiescent state, 
then studies of the blazar contribution to the EGB based 
on observed flux distributions such as that presented in 
Abdo et al. (2010c) could underestim,ate the number of 
quiescent blazars resulting in an underestimation of the 
blazar contribution to the EGB. Thus, the apparent dis- 
crepancy of the predictions of the two analyses could be 
the result of the variability of gamma-ray emission in 
blazars, the study of which could have implications for 
gamma-ray emission in blazars. 

The comparison of the two types of analyses could also 
provide insight in the cosmological properties of blazars. 
At 100 MeV, the Fermi- LAT angular resolution is ~ 5° 
(Atwood et al. 2009); hence, for blazar GLFs that predict 
large number densities, many blazars, particularly those 
on the faint-end of the source count distributions, that 
are. in principle, observable would not be distinguishable 
from other blazars (and thus are unresolved) and would 
not be included in source count distributions. In such 
a scenario, the faint-end of the blazar flux distribution 
would underestimate the number of blazars in a given 
flux bin, and the faint-end slope might appear flatter 


than it should be 9 , resulting in an underestimation of 
the blazar contribution to the EGB. Thus, the apparent 
discrepancy of the results of the two analyses could be 
the result of a large number density of blazars. 

In conclusion, we have demonstrated that even with 
the inclusion of spectral breaks, the collective spectrum 
of unresolved blazars reproduces well the spectrum of the 
Fermi EGB for . several break models. As such, we find 
that the possibility that the collective intensity of un- 
resolved blazars dominates the EGB is not excluded on 
spectral grounds, even if spectral breaks are shown to be 
common among Fermi blazars. Given the remaining con- 
troversy concerning the blazar GLFs, we conclude that 
it is, as yet, premature to rule out blazar dominance of 
the EGB. Furthermore, we have shown that relative con- 
tributions of the sub- populations of blazars required to 
fit the EGB spectrum are sensitive to the nature of their 
breaks; hence, the methodology we present in this paper 
can be used to constrain the GLFs of blazars. As models 
for the GLFs of blazars (accounting for the uncertain- 
ties outlined in this paper) and more data on spectral 
breaks become available, the study of the spectral shape 
of the blazar contribution to the EGB in light of the 
breaks can provide insight into the high-energy physics 
of blazars and their cosmological properties. 
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9 Another symptom (hut not proof) of the effect of source confu- 
sion would be a similar Fermi measurement of the EGB as that of 
EGRET at lower energies since at these energies, the angular reso- 
lution of Fermi-LAT is comparable to that of EGRET (Thompson 
et al. 1993; Atwood et al. 2009: Stecker & Venters 2011). Intrigu- 
ingly, the first few data points of the Fermi EGB do appear to 
be consistent with the Strong et al. (2004; determination of the 
EGRET EGB. 
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